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ABSTRACT 

o ■ 

o : 

I show here that the time-delay between the gamma ray burst proper and the 
onset of the afterglow is a sensitive function of the surrounding density, going from 
~ 10 s for typical disk densities to several hours for IGM densities. Since bulk Lorenz 
j^>| ' factor, explosion energy and environment density can be simultaneously determined 

^f-} ■ from observations of the time-delay and of the afterglow, afterglow observations alone 

may establish which fraction, if any, of all bursts resides inside galaxies, in their outer 
haloes, or outright in the IGM medium, subjecting the NS/NS merger hypothesis to a 
. definitive test. If bursts are due to collapse of massive stars, no matter how short the 

^ burst, the onset of the afterglow should be immediate. 

m ; 

, Subject headings: gamma rays: bursts 

o • 
23: 

Q-f 1. Introduction 

6; 

£3 , There are currently three favored models for the formation of GRBs: mergers of neutron star 

binaries (Narayan, Paczynski and Piran 1992), collapsars (or hypernovae, Woosley 1993, Paczynski 
j> \ 1998), and SupraNovae (Vietri and Stella 1998, 1999), but a clearcut test to which each of these 

could be subjected is still missing. This is so because the details of the fireball model are to a large 
extent independent of the details of the energy deposition mechanism. The situation has also a 



paradoxical side: even if we managed to observe the event that deposits the energy, the enormous 
expected optical depths (exceeding 10 10 ) would prevent us from deriving any information on the 
nature of the source. Gravitational waves can probably provide a signature of at least some of 
the scenarios, but detection of these from these distant sources does not look like an immediate 
prospect. 

For these reasons, many authors (Bottcher et a/. 1999, Meszaros and Rees 1998) have proposed 
that interaction of the burst with surrounding material might provide detectable signatures of 
star-forming environments. However, the possible discovery of iron emission lines (Piro et al., 
1999, Yoshida et al, 1999) from two bursts has highlighted the weak side of this possibility, i.e., 
that any conclusion is subject to the usual uncertainties associated with modeling complex and 
incomplete sets of observations (Lazzati et al, 1999, Vietri et al., 1999). I show in this paper that 
a surprisingly simple and robust test exists, for both the neutron binary merger and the collapsar 
scenarios, resulting from a trivial, but so far overlooked property of fireballs. 
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2. Time delays 

I consider below spherically symmetric fireballs. Since I am interested in phases of the evolution 
when the ejecta are still highly relativistic, the corrections due to the sideways expansion of a beamed 
model are expected to be small. The test is provided by the delay of the onset of the afterglow from 
the end of the burst proper. The burst duration is due exclusively to the total lifetime of the source, 
a quantity which we cannot currently predict, but which, given the large luminosities involved, is 
observed to lie within a hundred seconds. The radius at which the energy release occurs, at the 
contrary, is completely fixed by the requirement that the source geometry does not filter out the 
millisecond-scale variability which is so often observed: we must have R s h ~ 10 13 cm (Ruderman 
1975, Sari and Piran 1997). At these distances, the relative kinetic motion of the dishomogeneous 
wind is dissipated. The afterglow proper begins instead at a distance which is fixed solely by the 
density of the surrounding medium, and the overall properties of the wind: 



E A 1 / 3 / lcm' 3 \ 1/3 /300\ 2/3 
10 53 erg J \ n J \ 7 / 



R ag = A x 10 16 cm ( T7 ^- zz ) (±^_) ( — ) (1) 



where n is the number density of surrounding matter, and 7 the bulk Lorenz factor. At this radius, 
the inertia of swept-up matter becomes comparable to that of the ejecta, and slowdown begins. 
From this point on, the X-ray luminosity begins decreasing as t~ a , with a = 0.7 — 2.2, the complete 
range observed so far. The beginning of this easily recognizable phase occurs a time td after the 
end of the burst, with 

Rag ~ Rsh ^Ra l = ,, ( E \ ^ ( ± Cm~ 3 \ 1/3 /300 \ 8/3 



t d = ^ 2 "™ * = 15 si ^_ _ . (2) 
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The time-delay is a function of the product nj 8 , so that one may fear that large variations 
of n may be masked or mimicked by much smaller variations in 7; this degeneracy will be lifted 
shortly, but it should be remarked that that the ejecta Lorenz factor is limited from below to 
7 > 30 (Meszaros, Laguna and Rees 1993), and it is unlikely to exceed 7 ~ 300 by much. At 
the same time, considering for the moment 7 fixed, we see that the time delay td is a moderately 
sensitive function of the surrounding density: for a density typical of the interstellar medium in 
a galactic disk, n w 1 cm -3 , clearly td ~ 10 s. For long bursts, and in particular for the whole 
subclass of bursts accessible to BeppoSAX, this time-delay is comparable with the burst duration 
tb, and we would thus expect a superposition of the last parts of the burst and the beginning 
of the afterglow. However, the time-delay acquires interesting and measurable values when the 
burst takes place outside a disk: in a galactic halo, for which a baryonic representative density is 
around n w 10 -4 cm -3 , I find td ~ 5 minutes. In the case of the intergalactic medium, for which 
n « 10~ 8 cm' 3 , I find i d « Ah. 

So we see that the expected range of variation for 7 s , about 8 orders of magnitude, is similar 
to the expected range for n: if we find an independent way of fixing 7, then we may use the time- 
delay as a measure of whether the environment density really spans eight orders of magnitude, and 
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thus, in the end, of where the bursts dwell. That afterglow observations could fix the environment 
parameters had been noticed already in the literature; what is new here is the realization that the 
difference in environmental densities controls the time-delay as much as the ejecta Lorenz factor, 
and thus that the time-delay can be turned into a sensitive measure of the bursts' whereabouts. 

What happens between the end of the burst proper and the onset of the afterglow? Electrons 
accelerated at internal shocks, both by quasi-thermal and Fermi processes, cool much faster than the 
hydrodynamical expansion time-scale, (Sari, Narayan and Piran 1996). However, some electrons 
will be accelerated at the forward shock which separates the surrounding medium from shocked 
matter already raked up by the ejecta. Since the rate at which matter is raked up is roughly 
constant, we can take as a first approximation that the total emitted power is a constant, for which 
an upper limit can be found as follows. The shock slows down when the total swept up mass 
(times 7) equals the ejecta mass. Within this time interval, then, about half of the total energy is 
transfered to the swept up mass; if this were highly radiative, then the total emitted power would 
be of order E rj j- d which is much smaller (by the factor t^/^td, where £ « 0.1 is the ratio of relative 
to bulk kinetic energy of the wind) than the burst luminosity. However, this is an upper limit, 
because, though electrons are highly radiative, this estimate implies that they are able to radiate 
away also all of the thermal energy of the protons. Since it is generally assumed (and it can be 
checked in the observations of the afterglow of GRB 970508, Frail, Waxman and Kulkarni 2000) 
that only a fraction e eq « 0.1 of the protons' energy can be transferred from the protons to the 
electrons, this assumption is excessive by the same factor, e eq . The exact luminosity radiated during 
this coasting phase is difficult to predict exactly, since it depends critically upon rather uncertain 
plasma factors, like the efficiency e eq with which energy is tranferred from protons to electrons, the 
efficiency e# with which a near-equipartition magnetic field is built up, and the relative efficiency 
of quasi-thermal to Fermi processes. However, for the present aim, this exact value is unnecessary: 
suffice it to say that a reasonable expectation is that 

L<^- (3) 

What happens after shock deceleration has been described, for different values of the parame- 
ters, by Meszaros, Laguna and Rees (1993). Basically, a reverse shock will propagate at moderate 
Lorenz factors backward into the ejecta, dissipating an appreciable fraction of the kinetic energy 
of the shell within a timescale R ag /^ 2 c, which is the apparent time in which the shock crosses the 
ejecta shell, ~ td- Thus the maximum expected power is exactly a fraction t^/^td below that of 
the burst. Notice that, in the previous section, this was only an upper limit. This internal shock 
provides a new, second phase of emission, of total duration comparable to the light crossing time 
of the ejecta shell, td- 

After this peak luminosity we then expect the usual afterglow (oc t~ a ) to begin. A qualitative 
description of this can be seen in Fig. 1. During the afterglow, emission is due to synchrotron 
processes (Meszaros and Rees 1997), as evidenced by the simulataneous spectrum of GRB 970508 
that fits theoretical expectations so nicely (Galama et ai, 1998) and by detection of polarization 
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in the optcial spectrum of GRB 990510 (Covino et al, 1999, Wijers et al, 1999). In this case, it is 
well-known that the synchrotron peak frequency u m is given by (Waxman 1997a) 



rn 



where £53 is the explosion energy in units of 10 53 erg, td ay is the time of observations in units of 
days after the beginning of the afterglow, and 7 is the initial bulk Lorentz factor before the inertia 
of the swept-up matter begins slowdown. Using Eq. 2, it is convenient to rewrite this as a function 
of the time-delay t^: 



, m = 3 x 10 15 Hz (^) 1/2 (e eq /0.2f(e B /0.l)^E 53 t^ 2 (l^)^ . (5) 



As an example, for the time-delays discussed above, ten minutes after the beginning of the afterglow 
we expect the peak frequency to be hv m = 22 keV for the disk ISM, 250 eV for the halo, and 1 eV 
for the IGM. 

On the other hand, the intensity at v m is independent of the bulk Lorenz factor (Waxman 
1997b): 

but is a reasonably steep function of the environmental density. Eq. 2,5, 6 allow the simultaneous 
solution for the three unknown parameters E, n and 7, which are the only ones potentially varying 
by several orders of magnitude. In this way, we can lift the degeneracy present in the factor nj s . 
Basically, given a time-delay, a time tday after the new peak we may locate the afterglow peak at 
a frequency given by Eq. 5, and then the afterglow intensity (Eq. 6) will establish whether the 
time-delay is given a much smaller value of n or to a marginally smaller value of 7. In bands like 
the X-ray, if the time-delay is due to n <C 1 cm -3 , we basically do not expect any detectable flux, 
which might however be detectable in the optical/UV, while the afterglow is certainly detectable 
in the optical if a smaller value of 7 accounts for the time-delay. For this reason, experiments like 
SWIFT, which will provide simultaneous coverage in both optical/UV and in the X-ray, are of 
fundamental importance to interpret time-delays. 

Another interesting application of Eq. 1 occurs for the opposite case, in which the sur- 
rounding material is that of a pre-existing wind from the progenitor system. In this case, using 
n = M/4irr 2 m p v 00 , I find 

E Voo _ 14 E Voo W-SMpyr- 1 (300\ 2 

ag M 7 2 c 2 Cm i0 52 er<? 100 km M \ 1 ) ' l j 

From this it can be seen that only exceedingly short bursts, i.e., those lasting less than R ag /^ 2 c w 
0.03 s, can display a break in the X-ray emission: all other bursts (i.e., the near totality) must 
display an afterglow emission superimposed on the burst itself, and no break of continuity. 
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3. 



Discussion 



The relevance of the above arguments to the problem of the identification of the bursts' pro- 
genitors is clear. Bloom, Sigurdsson and Pols (1999) have shown that there is a well-defined 
distribution of expected distances of GRBs from their parent galaxies, if they originate from the 
merger of neutron stars. In particular, they showed that ~ 15% of all burst progenitors may be 
able to escape from their parent galaxy altogether, and that about 50% of all bursts will be located 
more than 8 kpc from their site of origin, nearly in the galactic halo. This feature is potentially 
verifiable also with optical investigations, but we know that even the most powerful ground-based 
telescopes are often inadequate to find the host galaxy, while HST searches, which up to now have 
scored a remarkable 5/6 success rate (Fruchter 2000), are slow, cumbersome, and, for some of the 
most distant bursts, cannot go deep enough in the luminosity function. The test I propose is instead 
simple, and a statistically meaningful database will be secured in the near future, since several space 
missions (AXAF, Hete II, Integral, XMM, SWIFT) will observe hundreds of bursts for several hours 
after the burst. Especially important among them is SWIFT, which will be responsible for most 
of these observations, and will provide continuous coverage of each detected bursts in both X-ray 
and UV/optical wavebands. From this database of perhaps 300 bursts observed continuously over 
several hours after the burst, we may expect ~ 50 bursts located in the IGM medium, and another 
~ 100 located in the haloes of their parent galaxies. Thus, within this scenario we expect about a 
third of all bursts to have delays of several minutes, and a few tens to have delays of a few hours. 
The flux should then peak again to a factor « tbl&d below that of the burst proper. 

The very same test can also be used a converso to check whether bursts are surrounded by the 
massive winds expected in all scenarios involving massive stars. In particular, it is expected that 
no time-delay will ever be observed, except perhaps for the shortest bursts, % <C 0.03 s. 

One may wonder what happens when only incomplete data are available, i.e., when the burst 
disappears shortly after tf, and is re-observed when already in the afterglow regime. This is the case 
of BeppoSAX which, because of instrumentational limitations, is incapable of providing significant 
information on the key silent period described here. It should be notice that the above estimates for 
the luminosity (L& at %, and L^td/^td at 2td) fall precisely on a t~ l slope. So, if further observations 
of the afterglow find an afterglow slope close to —1, extrapolation of the afterglow luminosity to early 
time will match closely with the burst luminosity, regardless of the existence of the silent period. 
Thus, when only incomplete data are available, only bursts with slopes significantly different from 
— 1 can be helpful in deciding whether a silent period exists. 

A corollary of the above tests comes from the observations of flares in the radio lightcurves. 
In fact, these flares are expected only from afterglow sources with radii smaller than R c (Goodman 
1997), where 



Here viq is the observing frequency in units of 10 GHz, d sc ^ pc the distance of the Galactic scattering 
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medium, assumed to be a simple screen, in units of 1 kpc, and SM is the Galactic scattering 
measure, scaled in terms of a typical Galactic value. It can be seen from the above that R c <C R a g 
every time that the burst is not located inside a galactic disk. Thus, for the ~ 150 sources mentioned 
above, which are expected to show measurable time-delays, we expect a fortiori no radio flares. 
Radio flares are instead allowed if bursts originate from collapsars. 

In summary, I have shown here that the X-ray luminosity has a silent period between the 
burst and the power-law afterglow, of duration 



15 s ISM 

5 min galactic halo (9) 
Ah IGM 



after which we expect a resurgence in the bolometric flux to a level a factor ~ tb/^td below that 
of the burst proper, from which the easily recognizable, power-law afterglow begins. Afterglow 
intensity at peak frequency depends only on the environmental density and total explosion energy, 
while the peak frequency depends on total energy, n and 7, allowing a simultaneous solution which 
fixes n (Eqs. 2, 5, 6). 
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Fig. 1. — Qualitative X-ray lightcurve of burst with measurable time-delay between burst proper, 
and afterglow. During % = 10 s, there is the X-ray- bright burst; during the next t c < the 
electrons accelerated at the internal shocks cool, and the X-ray luminosity decays to a small value; 
for t < td = 10 3 s, the X-ray flux is due to the small fraction of the (small) radiated luminsoity 
which ends in the X-ray band: most of the flux will be in the 7-ray band. When t « td = 10 3 s, the 
forward shock begins decelerating, driving a reverse shock into the ejecta; this causes a resurgence 
of the X-ray luminosity on a time-scale ~ td, up to a flux level a factor below that of the 

burst. At a time ~ td, the afterglow proper begins, first driving a reverse shock into the ejecta 
(which is responsible for the new peak), then beginning the the power-law-like decay. Notice that 
the normalization of the afterglow (Eq. 5, 6) is unrelated to that of the peak. 



